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PACS 73 . 23 . -b - Electronic transport in mesoscopic systems 
PACS 73 . 23 . Ra - Persistent currents 

PACS 42. 65. Re - Ultrafast processes; optical pulse generation and pulse compression 

Abstract. - We study theoretically the charge polarization and the charge current dynamics of a 
mesoscopic ring driven by short asymmetric electromagnetic pulses and threaded by an external 
static magnetic flux. It is shown that the pulse-induced charge polarization and the associated 
light- emission is controllable by tuning the external magnetic flux. Applying two mutually per- 
pendicular pulses triggers a charge current in the ring. The interplay between this nonequilibrium 
and the persistent currents is investigated and the conditions under which the pulses stop the 
persistent current are identified. 
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INTRODUCTION. - Quantum structures with a 
ring confining geometry have served over the years as a 
paradigm for the demonstration of quantum-interference 
phenomena such as the Aharonov-Bohm effect and the 
persistent (equilibrium) current [1-4] that emerges if the 
ring is pierced by a static magnetic flux. The persistent 
current has also been considered when a time-dependent 
magnetic field with a static component [5-12] being ap- 
plied, and it has been found that the current diminishes 
when the static component vanishes. Another possibility 
for generating (nonequilibrium) current is to irradiate the 
ring with circular polarized light, [13-15] and references 
therein. In principle, this approach is feasible for a size- 
quantized system and if the radiation frequency is tuned 
such that the rotating-wave approximation is applicable 
and counter-rotating contributions to the current are thus 
negligible. Hence, this approach is expected to be par- 
ticularly useful for molecular ring structure, cf. ref. [15] 
and references therein. Another approach for current gen- 
eration which does not rely on the level quantization and 
resonant excitations, is based on an asymmetry of the elec- 
tric field amplitude of the applied pulse [16,17]. As exper- 
imentally demonstrated [18-21] the optical cycle of such 
asymmetric pulses consists of a short and strong half cycle 



followed by a much longer and weak half cycle of an oppo- 
site polarity. The charge dynamics is mainly triggered by 
the strong , short half cycle [16,17,22,23] and hence these 
pulses are referred to as half-cycle pulses (HCP). In gen- 
eral, HCPs are capable of inducing nonequilibrium electric 
dipole and direct charge currents [17,22,23], even if the 
system is not lacking an inversion symmetry, in contrast 
to the well-known photovoltaic effect [24,25]. For isolated 
rings (and without an external magnetic flux), applying 
one linearly polarized HCP induces a charge polarization. 
Applying in a perpendicular direction to the first pulse a 
second HCP leads to a charge current that can be tuned by 
varying the properties and the time-delay between the two 
pulses. The question which has not yet been addressed and 
will be of concern here is how the HCPs-induced charge 
polarization and the charge current are influenced by the 
persistent current generated by a static magnetic flux. As 
shown below, the HCP-induced nonequilibrium charge po- 
larization and the associated emitted radiation are tun- 
able by changing the magnitude of the static magnetic 
flux. The nonequilibrium current triggered by two time- 
delayed HCPs can be tuned to cancel the magnetic-flux- 
induced equilibrium current offering thus a possibility for 
time switching of the persistent current on a picosecond 
time-scale. 
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Theoretical model. — For the sake of simplicity, we 
consider here an isolated single channel ID ballistic ring 
with N electrons and radius po at low temperatures. This 
treatment is appropriate if the width and the height of 
the ring are smaller than both the ring radius and the 
Fermi wavelength of the carriers. The generalization to 
the multi-channel case can be done along the lines of 
Rcfs. [16,26] where it is shown that the multi-channel case 
does not change qualitatively the predictions of the ID 
model. The ring is pierced by an external magnetic flux 
(j> and is subjected to a sequence of short, linearly polar- 
ized half-cycle pulses. The light propagation direction is 
perpendicular to the plane of the ring. The Hamiltonian 
describing the dynamics of the system is given by 



H(t)=H + V int , 



(1) 



where 



H 



d_ 



. r/ ,, l + 5 V int = eE(t) ■ P- (2) 



Here e is the elementary charge, E(i) stands for the time- 
dependent electric field of the applied pulses and the vec- 
tor p — (p n cos9,posin9) describes the electron position. 
m* and <f) n are the electron effective mass and the flux 
quantum, respectively. In what follows we expand the 
field operator ^>(9,t) on a basis of eigenstates <p n (0) of 
H , i.e., 

*(M)=^c„(t) fc («), H oVn (9) = E n °^ n (9). (3) 



The components of the time-dependent density matrix are 
expressed as p mn {t) = (cln(t)c n (t)}. Here (...) denotes the 
expectation value which is taken with respect to the initial 
state of the system. 

Now we turn to the specification of the employed light 
pulses. The duration Td of the experimentally available 
HCPs is in the subpicosecond time scale. On the other 
hand, for typical ballistic mesoscopic rings (po ~ lpm) 
the time r F a particle at the Fermi level needs for a round 
trip is tens of picoseconds, meaning that r<j <C t f . In 
this case, the so-called impulsive approximation (IA) can 
be employed to describe accurately the dynamics of the 
system [27]. Within the IA the action of the HCPs is en- 
compassed in the matching conditions (4) for the density 
matrix at the instances before and after the pulses are ap- 
plied while the system propagates in a pulse-free manner 
at any other time [28] (note that in absence of the pulses 
the magnetic flux is still present). To be specific, let us 
consider the case where two HCPs that are linearly po- 
larized in the x and y directions are applied respectively 
at the times t\ = and t 2 = r. For the density ma- 
trix we find within the IA that (we use the abbreviation 



Wn{9)) = |n» 
/W(0+) = ^(n| e - MlCOse |0(fc|e MlCOse |m) Pfei (0-); 



k,l 



Pmn(r + ) = ]>>|e- M2Sin0 |O(% M2Sin VWO- 



(4) 



Here the dimensionless quantity ctj characterizes the ac- 
tion transferred to the system 



ttj = pop, IK 
while pj describes the strength of the pulse Ej 



Pi 



"J 

Jt. 



tj+T d 



Ej(t)dt. 



(5) 



(6) 



The subindcx j = 1,2 refers to the first and second pulse 
(with the same duration Td). We remark here that for a 
different propagation direction of HCPs, e.g. in the plane 
of the ring, and/or for the strong excitation regime (a 3> 
1) the time-dependent magnetic-field component of the 
pulse may in general affect the dynamics of the electron 
system and should be included in the Hamiltonian (1). 
This is to be contrasted with the case discussed in this 
paper where the magnetic field associated with the HCP is 
moderate and lies in the plane of the ring (where electrons 
are confined). 

Within the relaxation time approximation the equation 
of motion for the density matrix reads 



do E (0) - E {0) 

u f J mri . -*^m 



3mn r /~_\-| ' $mn 

rj-y^ [Pmn ~ Pn7i(v J J 7j~, Pmn-> 



(7) 



where 8 mn is the Kronecker symbol and T\ and Ti are the 
relaxation and dephasing times, respectively. Eqs.(4)-(7) 
together with the initial condition that before the applica- 
tion of the pulses the system is in a thermal equilibrium, 
i.e., 



Pmn{0 ) = 25mn 1 + e k B 



(8) 



determine completely the time evolution of the density 
matrix. The factor 2 in the Fermi-Dirac distribution (8) 
accounts for the two-fold spin degeneracy and the chemical 
potential rj is set by the requirement of the conservation 
of the particles' number N . 

We are particularly interested in the study of the charge 
polarization and the currents induced in the ring. The 
charge polarization along the x axis is determined by the 
^-component of the dipole moment operator. After some 
algebra one obtains for the charge polarization 



p(t) = — ep Rc 



E. 



(9) 
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This result is remarkable. It tells us that the time de- 
pendence of the charge polarization does not depend on 
the time evolution of the diagonal components of the den- 
sity matrix, i.e., by tracing the evolution of the induced 
polarization one can obtain physical information that is 
exclusively due to dephasing. 

Using the first of eqs. (4) and eqs. (7)- (9) and perform- 
ing the same steps as in ref. [22] we can write the charge 
polarization of the ring after the application of the first 
HCP as 



-epoai©(*)[Jb(fi) + MO)] sin 

~47rf 



2vrf 



-t/T 2 



^Pnn(0 )C0S 



(10) 



where 0(x) denotes the Heaviside step function and Ji(x) 
is a Bessel function of the order I. Furthermore, we intro- 
duced the quantities 



Q = ai ^2 - 2cos(47rf/i p ), 

t p = 47TTO* Pq/K. 



(11) 

(12) 

The charge current density operator can be expressed in 
terms of the field operators (3) and the angular component 
of the vector potential Ag = (f>/(2np ) as follows 



3 = 



ieh 
2m*pl 



pomrc 



86 
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(13) 



Upon mathematical manipulations we find for the charge 
current 

I(t) = «j» = IoJ2( n + MMPnnit). (14) 



Here 



I = -eh/(m*p ) 



(15) 

sets the scale of the current magnitude. The double 
bracket ((•••)) stands for angular integration and expec- 
tation value computation. In contrast to the charge polar- 
ization, the charge current depends only on the diagonal 
elements of the density matrix. Therefore, dephasing and 
relaxation can be independently investigated by tracing 
the charge polarization and the charge current, respec- 
tively. 

Exploiting Eqs. (4)-(7) and (14) we deduce after some 
algebra the following relation for the total charge current 1 



where 



I(t) = Ip ers +Q{t- T)I dyn (t), 
Ipers ^0 



(16) 

(17) 



1 Self-inductive effects turned out to be completely negligible in 
comparison to Ipers and Idyn(t)- 



is a static (persistent current) component induced by the 
magnetic flux and 



h V n{t) = ^%(r)e-^ 
epo 



(18) 



is a dynamical component due to the action of the pulses 
and is determined by the polarization fj,(r) of the ring at 
the time the second pulse is applied [16]. Note also that 
the dynamical component of the current only appears af- 
ter the second pulse is applied, i.e., it is necessary to ap- 
ply two orthogonal linearly polarized HCPs in order to 
produce an additional clock-wise anti-clock-wise symme- 
try breaking [16] (in case of the persistent current this 
symmetry break is brought about by the magnetic-flux 
induced time-reversal symmetry breaking). 

It follows from eq. (16) that the static part of the total 
charge current induced in the ring does not depend on the 
application of the HCPs and is not modified by them (i.e., 
the persistent current is really persistent). On the other 
hand, the dynamical component Id yn depends on both the 
external magnetic flux and the pulse parameters because, 
as shown below explicitly, the dipole moment p,(r) is a 
function of <f>. 

An important consequence of eqs. (16) and (18) is that 
a magnetic-flux induced persistent current in a ballistic 
ring can be temporally stopped by applying HCPs, for all 
the parameters determining the sign and the magnitude of 
Idyn are known [cf. eqs. (5), (15) and (19)]. We note that 
the build-up time of Id yn (which determines the switch- 
off time of I, i.e. the cancellation of I per s) is set by the 
duration of the HCP which can be in the subpicosecond 
regime. As Id yn relaxes according to eq. (18) I pe rs emerges 
again, however, upon the application of another sequence 
of HCPs we can switch off again the total current (16). 

It should be remarked here that the decomposition of 
the charge current in a statical and a dynamical compo- 
nent is valid within the relaxation time approximation em- 
ployed here. It might well be that the inclusion of memory 
and non-linear effects prevents such a simple structure of 
the current. A definitive conclusion on this point would 
require the treatment of two-particle and higher order cor- 
relations that is a challenging task for inhomogeneous dis- 
tributions (non-diagonal density matrix elements). On the 
other hand, for weak excitations and small temperatures 
we expect nonlinear effects (in the field strength and in 
relaxation) to be marginal and the correlation effects on 
the relaxation to be of less importance (due to the same 
reasons as for the ground state). 

Explicit results and numerical demonstrations. 

For zero temperature we were able to perform the sum in 
the eq. (10) analytically to obtain 

p(t) = -epoai6(t)[J (0) + J 2 (Q)]s(i) e -'/ T2 , (19) 
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Fig. 1: The time dependence of the charge polarization p(t) as a function of the magnetic flux cf> for various values of the 
pulse strengths that result in different values of the transferred action ai (cf.Eq.(5)). For the sin-square HCP having duration 
Td — 1 ps and the ring radius of ro = 1 A™ the value of the transferred action ai = 0.1 («i = 1, a = 10) corresponds to the peak 
electric field strength Eq = 1.32 V/cm (Eo = 13.2 V/cm, Eo = 132 V/cm). The time is expressed in units of t p = 4irm* Po/h. 
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These results apply for 0/0o £ [—1/2,1/2], otherwise 
s(t; 0/0o) is given by the periodicity condition s(t; 0/0o + 
l)~s(t; 0/0o). In the limiting case = eqs. (20) simplify 
to the results obtained in ref. [22]. Inserting the expres- 
sions for ^{t) given by eqs. (19) and (20) at t = r into 
eq. (18) yields an analytical expression for the dynamical 
part of the current. The persistent part of the current was 



calculated previously in ref. [29] . 

For an illustration we performed explicit calculations 
for the case of a ID GaAs ring with N = 100, radius 
pa = 1 fim at zero temperature. The time dependence 
(ignoring the effects of dephasing [30]) of the induced 
charge polarization is shown in fig. 1 for different val- 
ues of the external magnetic flux and pulse strengths 
yielding different ol\ [cf. eq. (5)]. The time is expressed 
in units of the characteristic time t p . For certainty we 
assume here that the time profile of the HCP electric 
field is given by E(t) = E sin 2 [irt / T d ] for t e (0, Td) and 
E(t) = outside of this time interval. In this case we have 
oi\ — poeE Td/(2h). As seen from fig. 1 different patterns 
of the charge polarization evolution can be tailored by 
changing the applied flux and/or the pulse strength. This 
observation is of particular relevance since the emission 
properties of the ring are determined by the time oscilla- 
tions of the charge polarization [17,31]. Thus, the system 
proposed here could serve as a source of electromagnetic 
radiation with magnetic- flux controllable properties. An- 
other remarkable fact revealed by fig. 1 is the presence 
of beating behaviour of the charge polarization dynamics 
in the low excitation regime as well as collapses and re- 
vivals of the polarization at higher excitations signifying 
the important role of quantum interferences [17]. 
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Fig. 2: Dependence of the peak current / on the magnetic flux 
<f> and the time delay r between the pulses for different values 
of the pulse strengths: a) ai = «2 = 2, b) ai — cti — 5 and c) 
Ql = Q2 = 10. 

Now we consider the current dynamics. As an example, 
fig. 2 shows the dependence of the peak current on the 
magnetic flux and the time delay between the two pulses 
that trigger Id yn - The peak current is given in units of Io 
(e.g., for a ring with po = 1 ^m and N = 100 we have 
Io = 3 nA). The total charge current has an oscillatory 
behavior when scanning cj> and r. In particular, both the 
magnitude and the sign of the peak current can be changed 



by appropriately varying (f> and r. As we outlined above, 
there exist regions in which the total induced current van- 
ishes due cancellation of the persistent and the dynamical 
components of the current. As it is seen from the figure 
the region where the total current is switched off is con- 
trollable by sweeping <f>, r and/or the pulse strengths. All 
of these parameters are externally tunable and experimen- 
tally feasible with current technology. 

Strictly speaking, the derived analytical formulas, (16)- 
(18) and ((19-20)), are only justified for the weak excita- 
tion regime. In the strong excitation regime the nonlinear 
terms in relaxation can couple diagonal and non-diagonal 
components of the density matrix, and in this way the dy- 
namics of the polarization and current. However, if the 
relaxation is slow (and we can expect it to be so for low 
temperatures) it hardly affects the results presented in 
the numerical illustrations also for the strong excitation 
regime (a = 5, 10). 

Finally we note that the circulating current results in an 
induced (orbital) magnetization (M(t) ~ icp^I {£)"). The 
current temporal controllability is reflected in a picosecond 
switching of the ring magnetization. The above finding 
can be generalized, particularly in view of potential appli- 
cations, to the case of trains of HCPs pulses and/or arrays 
of mesoscopic rings [16,17]. The first case may serve as 
a sophisticated source of ultrashort magnetic pulses while 
the second allows the design of artificially microstructured 
materials with magnetic properties locally controllable. In 
view of applying the present scheme to nano and molec- 
ular rings [15] we note that this regime is accessed by 
changing appropriately the parameters of HCPs: The ap- 
plicability of the above theory is based on the assumption 
Td -C t p which implies the use of femtosecond HCP for 
nano-size rings. Such pulses and even attosecond HCPs 
are currently under discussions [32]. 

Conclusions. — In summary, we showed that the 
time-dependent charge polarization and charge currents 
can be generated in mesoscopic rings threaded by a mag- 
netic flux and subjected to a sequence of asymmetric 
electromagnetic pulses. The time dependence of the in- 
duced polarization shows different patterns that can be 
controlled by appropriately adjusting the magnetic flux 
and the pulse parameters. Quantum interferences result in 
beatings and revivals of the charge polarization dynamics. 
The total charge current induced in the ring can be de- 
composed into a static and a dynamical components. The 
static component is associated with the persistent current 
and does not depend on the parameters of the pulses. The 
dynamical component is expressible analytically and can 
be tuned to cancel the persistent current allowing thus a 
picosecond switching of the total current in the ring. 
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